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Abstract 



This paper presents the exact expressions of the transition probabilities of some 
non-determinantal Bethe ansatz solvable interacting particle systems: the two- 
sided PushASEP, the asymmetric avalanche process and the asymmetric zero 
range process. The time-integrated currents of the asymmetric avalanche pro- 
cess and the asymmetric zero range process are immediate from the results of 
the asymmetric simple exclusion process. 



1 Introduction 



The Bethe ansatz technique is one of the methods widely used to find exact solutions 
of interacting particle systems. In physics the Bethe ansatz technique was origi- 
nally used to obtain the eigenvalues and the eigenvectors of the Hamiltonian of the 
one- dimensional quantum spin- 1/2 chain [22J. However, from the fact that Markov 
matrices of some interacting particle systems are a similarity transformation of the 
Hamiltonians of quantum spin chains one can study dynamical properties of the in- 
teracting particle systems by using the Bethe ansatz [5]. 
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Among the Bethe ansatz solvable models, the asymmetric simple exclusion pro- 
cess (ASEP) is a paradigmatic model in non-equilibrium statistical mechanics. A 
sub-model of the ASEP, the totally asymmetric simple exclusion process (TASEP) 
has been connected to random matrix theory due to the seminal works [3j [11] (for 
review, see [7J [13]) and some results on the TASEP have recently been extended to 
the ASEP [24-28]. In this direction, the mapping between the TASEP and the cor- 
ner growth model based on the RSK algorithm and the Bethe ansatz technique have 
been mainly used. Schiitz [20] found the transition probability of the TASEP with N 
particles and expressed it as an N x N determinant whose entries are hypergeomet- 
ric functions. This exact solution of the Kolmogorov forward equation is a starting 
point of time-integrated currents of the system and related works. Indeed, the current 
distribution of the TASEP and its large time asymptotics were well studied by us- 
ing the determinant representation and some properties of entries of the determinant 
[5j [TUl HH HZl HE] ■ This determinantal structure of the TASEP is inherited from the 
S'-matrix that describes a two-particle interaction in the TASEP. The S'-matrix in 
the ASEP is 

s i + A&zk (£,,<;, ec). (i.i) 
v + qiaip - £ a 

In case of the TASEP, that is p — 1 and q = 0, (II -ID is separable in the sense that 
the numerator has only £^ variable and the denominator has only £ a variable. This 
separability of variables makes it possible to express the transition probability as a 
determinant. However, as one sees, the ASEP does not have this property and thus 
the determinantal structure is no longer in the ASEP. 

The determinantal structure appears in the drop-push model [TH], too. In the 
drop-push model a particle can jump to the right neighboring site with rate 1 even 
if it is already occupied, pushing all the right neighboring particles by one. That is, 
representing a particle by A and an empty site by 0, the process 

A---A0->0A---A (1.2) 



occurs with rate 1 for all n G N when the exponential clock of the leftmost particle 
in (11.21) rings. Even if we allow particles to jump in both directions, the resulting 
process is still determinantal if the pushing effect is one-sided. Indeed, Borodin 
and Ferrari [4J have shown that limiting processes are the Airy x process and Airy 2 
process for flat initial condition and step initial condition, respectively. For the ASEP, 
which is not determinantal, in a breakthrough [24-26] Tracy and Widom obtained 
the distribution of time-integrated current and its large-time asymptotics for step 
initial condition. Tracy and Widom found first the transition probability of the 
ASEP with N particles and then obtained an integral formula for the distribution 
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of time-integrated current by summing the transition probability over all possible 
configurations. A novel point of Tracy and Widom's transition probability is that it 
is given by an integral representation over arbitrarily small contours or large contours 
keeping the forms of S'-matrices in the integrand. Thus, one can find the probability 
distribution of the mth particle's position by using various combinatorial identities 
related to the S-matrix of the ASEP. 

In this paper we find the exact expression of transition probabilities of other 
Bethe ansatz solvable interacting particle systems which are not determinantal. The 
transition probabilities are a starting point of working on the currents of systems as 
we saw in Tracy and Widom's works. The models we will consider are the two-sided 
PushASEP, the asymmetric avalanche process (ASAP) and the asymmetric zero range 
process (AZRP) with constant rates. We shortly introduce the models. If the process 
( II ,2p occurs with rate r n dependent on n, the process is no longer determinantal, and 
in this case, r n is fixed to 

' (A + /i = l) (1.3) 



1 + X/ n + (A//i) 2 + ■ ■ • + (A/a*)» 



to permit the Bethe ansatz solvability [lj. This model which is called the generalized 
totally asymmetric exclusion process or the one-parameter family of asymmetric ex- 
clusion process was first suggested by Alimohammadi, Karimipour and Khorrami pp. 
The generalized totally asymmetric exclusion process interpolates between the TASEP 
(A — Y 1) and the drop-push model (A — > 0). Later this model was extended to the 
asymmetric case called the two-parametric family of asymmetric exclusion process [2]. 
In the asymmetric case the process is 

A_-_A0 ^ A---A (1.4) 

n n 

with rates pr n to the right and ql n to the left where r n is (11.31) and 

ln= i + n/\ + w\)2 + ... + {(1 /x)»-i> ( A + / i = 1 )- ( L5 ) 

When p, q ^ 0, this asymmetric model interpolates between the model with the push- 
ing dynamics on both sides (A,yU ^ 0) and the model [I] with the TASEP dynamics 
on one side and the pushing dynamics on the other side (A (or //) = 1). In this paper 
we will call this two-parametric family of asymmetric exclusion process simply the 
two-sided PushASEP and assume that A and [i are nonzero. 

In the ASAP [151 EE] particles jump to the right or to the left with rate p and q, 
respectively, after an exponential time. When particles try to jump to the site already 
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occupied, an avalanche starts immediately. The rule for the avalanche is as follow^: 
If n (> 2) particles are at site x during an infinitesimal time, then n particles jump to 
x + 1 all together with probability fi n (fi n ^ 1) or with probability A n , n — 1 particles 
jump to x + 1 all together and one particle stays at x where /i n + A„ = 1. For example, 
for a two-particle system with p — 1, q — 0, if particles are initially at x and x + 1 
and the exponential clock at x rings at time t, then the resulting configuration is one 
of (x + n, x + n + 1), n G N and the probability of finding two particles at x + n and 
x + n + 1 is A/i n_1 where /x = /12 and A = A2. Regarding the Bethe ansatz solvability, 
it is known [151 IB] that /i n must be given by 

1 _ (—u)ri-l 

At n = /i f-^ , n>2. (1.6) 

1 + // 

As a special case, if fi = —p/q, the process becomes the ordinary ASEP [IB] , and it 
is easily seen that the ASAP with = is equivalent to the two-sided PushASEP 
with = 1 and so the process is determinantal. In this paper we consider only the 
case < /i < 1. 

The asymmetric zero range process (AZRP) is also not determinantal (p, q 7^ 0). 
The dynamics of general AZRP is governed by the following laws; if a site x is occu- 
pied by n particles, one of the particles at x leaves the site x at the rate g(n). In other 
words, the jumping rates g(n) depend on the number of particles at the departure 
site. The particle that leaves x chooses a target site y with probability p(x,y). In 
this paper, we assume that p(x, x + 1) = p and p(x, x — 1) = q with p + q = 1 and 
g(n) = 1 for all n. 

The paper is organized as follows. In section 2, we discuss the forward equations 
and the ^-matrices of the models as a background. In section 3, we provide the 
main result, the transition probabilities (I3.28P of the models, and using the obtained 
transition probabilities, we revisit the map between configurations of the AZRP and 
configurations of the ASEP [6j [9l [12l [21] . Also, we shortly discuss the current distri- 
butions of the ASAP and the AZRP. 



2 Preliminaries 



In this section we review the standard method by the idea of Bethe to discuss the 
^-matrices of the models that characterize each model. The Bethe ansatz solvable 
models in this paper with a finite number of particles are continuous-time Markov 

1 It takes an infinitesimal time for the avalanche to end. 
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processes with countable state spaces. A configuration of the models is specified 
by each particle's position and we denote the configuration by X = (xi, • • • ,x^) € 
Z N . Let Py(X;t) = P Y ( x±, ■ ■ ■ ,xjsr',t) be the transition probability from the initial 
configuration Y to a configuration X at time t and we call the domain of Py (X; t) as 
a function of X the physical region of the model. The physical region of the two-sided 
PushASEP and the ASAP is 

{(xi, • ■ ■ , xn) £ Z N '■ x\ < ■ ■ ■ < xn} 

and the physical region of the AZRP is 

{(xi, • • • , x N ) G Z N : Xi < ■ ■ ■ < x N }. 

Depending on X, the Py(X;t) satisfies different forward equations. For example, in 
the TASEP with two particles, Py{xi,X2',t) with x\ < x<i — 1 and Py{xi,X2',t) with 
X\ = x 2 — 1 satisfy 

^- P Y (x 1 ,x 2 ;t) = P Y (xi - l,x 2 ;t) +P Y (x 1 ,x 2 - l',t) - 2P Y (x 1 , x 2 ; t) (2.7) 
at 

and 

— Py(xi, x 2 ; t) = Py(xi - 1, x 2 ; t) - Py(xx, x 2 ; t), (2.8) 

respectively. However, introducing a function u(xi, x 2 \ t) on Z 2 x [0, oo) that satisfies 
(12. 7p for noninteracting particles and a boundary condition at the boundary of the 
physical region, two differential equations (12 .7p and (12. 8p are combined as a single 
differential equation with the boundary condition. That is, for u(xi,x 2 ;t) on Z 2 x 
[0,oo), 

— u(xi,x 2 ; t) = u(xi - 1, x 2 ; t) + u(x x ,x 2 - 1; t) - 2u(xi, x 2 ; t) (2.9) 

with 

u(x, x; t) — u(x, x + 1; t) (2-10) 

is in the form of (12. 7p when X\ < x 2 — 1 and in the form of (12.81) when x — X\ — x 2 — 1. 
Hence, if u(xi,x 2 ;t) satisfies (I2.9p . (I2.10p and the initial condition 

u(xi,x 2 ;0) = 5 yi (x 1 )5y 2 (x 2 ) when x\ < x 2 , 

then u(xi, x 2 ; t) = Py(X; t) with X = (xi, x 2 ) and Y = (yi, y 2 ) in the physical region. 
As a Bethe ansatz solvable model, the TASEP with iV particles (N > 3) does not 
need new constraints except the one in the form of (I2.10p . Hence the solution of 

d N 

— u(X;t) = ^ (u(x u ■ ■ ■ ,Xi_i,Xj-l,x i+ i, • • ■ ,x N ;t)-Nu(x 1 , ■ ■ ■ ,x N ;t)j (2.11) 

i=l 
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with the boundary condition 

u(xi, ■ ■ ■ , x h Xi + 1, • - • , x N ; t) = u(xi, ■■ ■ , Xi, Xi, ■ ■ ■ , x N ; t) (2.12) 

and the initial condition 

u(X; 0) = Sy(X) when x\ < ■ • • < x^ 

is the transition probability from the initial configuration Y to a configuration X at 
time t. 

Similarly, the forward equation of the drop-push model with N particles is (12. lip 
and the boundary equation is 

u (xi, ■ ■ ■ , Xi — 1, Xi, ■ ■ ■ , xn] t) = u(xi, ■ ■ ■ , Xi, Xi, ■ ■ ■ , xn; t). (2.13) 

Alimohammadi et al. [H [2] combined f !2.12j) and fl 2 . 1 3 [) in the form of 

u(xi, ■■■ ,x u xi, ■ ■ ■ ,x N ;t) (2.14) 
= im(x u ■■■ ,Xi-l,Xi,--- , x N ; t) + Au(xi, ■ • • ,Xi,Xi + l,--- ,x N ;t) 

and have shown that (I2.14p describes the process (11.21) with rates (11.31) or the process 
(11.4j) in the asymmetric case with rates pr n and ql n given by (II. 3p and (II. 5p . It is 
possible to approach Alimohammadi et al.'s model in a different point of view. Let us 
assume that we have an interacting particle system of which process is (11.41) with rates 
pr n and ql n to the right and to the left, respectively. These rates are to be determined 
so that the model is Bethe ansatz solvable. Let us consider the two-particle system. 
Setting n = h — 1, the forward equation is for u{x\, x% t) on Z 2 x [0, oo) 

4- u(xi, x 2 ; t) = pu(xi - 1, x 2 ; t) + pu(xi,x 2 - 1; t) (2.15) 
at 

+qu(xi + 1, x 2 ; t) + qu(xi,x 2 + 1; t) — 2u(x±, x 2 ; t), 

and the boundary condition is 

pr\u{x, x; t) — pr 2 u{x — 1, x; t) — p{r\ — r 2 )u(x, x + 1) (2-16) 
= —ql\u(x + 1, x + 1; t) + ql 2 u(x + 1, x + 2; t) + q(li — l 2 )u(x, x + 1) 

for all x G Z. However, instead of (12.161) . we consider a sufficient condition for it. To 
do so, we set the right hand side and the left hand side of (I2.16P equal to zero, and 
letting X = l 2 /h, = r 2 /r 1 and A + /i = 1, then both equations reduce into one single 
equation 

u(x, x; t) = fiu(x — 1, x; t) + \u(x, x + 1; t). (2-17) 
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Hence, if u(xi,x 2 ',t) satisfies (I2.17p . it also satisfies (I2.16p . This is the boundary 
condition suggested in [TJ [2] as a combination of the boundary conditions of the 
TASEP and the drop-push model. Hence, we take (I2.14p as the boundary condition 
for a two-particle sector (when only two particles are adjacent) of an iV-particle 
system. It has been shown [1] that (12.141) implies the boundary conditions for a 
general iV-particle system if r n and l n are given by (11.31) and (ll.5p . and so the model 
is Bethe ansatz solvable. Applying (I2.17P to the Bethe ansatz solution 

A (XI (X2 I A (X2CX1 

^12^1 ?2 + A 2Kl ?2 

generates the 5-matrix 

s " a -~T a '^ + xu,-^-T a ' s ^ (2 ' 18) 



In the ASAP with two particles if x\ < x 2 — 1, then the forward equation is (I2.15p . 
and if X\ — x 2 — 1 = x, the forward equation is 

— u(x,x + l;t) = pu(x — 1, x + 1; t) + qu(x, x + 2; t) (2.19) 
at 

oo 

+X(p + qyi) /i n ~ 1 w(x — n, x — n + 1; t) 

n=l 

—u(x, x + 1; t) — (j> + qfi)u(x, x + 1; t). 

Thus the boundary condition 

oo 

pu(x, x] t) — p\ / u n_1, u(x — n, x — n + 1; t) (2.20) 

n=l 

oo 

= —qu{x + 1, x + 1; t) + qX fi n u(x — n, x — n + 1; t) 

n=l 

is obtained by subtracting (I2.19p from (12.151) with x\ = x,x 2 — x + 1. However, 
setting both sides of (I2.20p to be zero, (12.201) is satisfied if u(xi, x 2 ; t) satisfies 

oo 

u(x, x;t) — X fi n u(x — n — 1, x — n; t). (2-21) 

n=0 

It can be shown that (12.211) is equivalent to 

u(x, x\ t) = Xu(x — 1, x; t) + fm(x — 1, x — 1; t) (2.22) 
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by applying ( I2.22p to itself recursively. Hence we take 

u(x u ■■■ ,Xi,Xi,--- , x N ; t) (2.23) 
= Xu(xi, ■ ■ ■ , Xi - 1, Xi, ■ ■ ■ , x N ; t) + fiu(xi, ■■ ■ ,Xi - l,Xi - 1, - ■ ■ ,x N ;t) 

as the boundary condition for a two-particle sector of the ASAP with N particles, 
and the resulting S-matrix is 

Spa = -» + \£ a -te a (2 - 24) 

The boundary condition (I2.23P implies the boundary conditions for iV-particle inter- 
actions when /i n is given by ( II. 6p [T5| [16] . 

In case of the AZRP with the assumptions mentioned in section 1, the S'-matrix 
generated by the boundary condition 

u(x 1 , ■■■ ,Xi,Xi,--- ,x N ) (2.25) 
= pu(xi, ■■ ■ ,Xi,Xi - 1, • • • ,x N ) + qu(xi, ■■ ■ ,Xi + l,Xi, ■ ■ ■ ,x N ) 

is 

o _ v + gjajp - _ ia q t 

Setting up the problem with N particles, for all three models we have the forward 
equation 

N 

u(X;t) = ^ (pujxi, ■■■ , Xj-i, Xi ~ • • • ,x N ;t) (2.26) 



d 
dt 



i=i 

+ qu(x!, ■■■ , Xi-x, Xi + 1, x i+ i, ■■ ■ , x N ; t) - Nu(X; t) 
and the initial condition 

u(X;0) = 5 Y (X) (2.27) 

which holds in the physical region of the corresponding model, and finally, the bound- 
ary conditions ([2TTI]) . ( |2T23|) and (12725] for the two-sided PushASEP, the ASAP and 
the AZRP, respectively. 



3 Transition Probabilities for iV-particle system 

In case of the two-sided PushASEP, the transition probability was expressed [2] as 
an integral over the unit circle, and thus S'-matrices in the integrand are not exactly 
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in the form of (12.181) to satisfy the initial condition. However, as a first step to 
compute time-integrated currents we need the integral representation over a contour 
with sufficiently small radius or over a contour with sufficiently large radius so that 
a certain geometric series converges [23]. So, in this section we find the transition 
probabilities in the form of the one of the ASEP in [21]. 

As usual, the energy in the Bethe ansatz solution is £tv(£) = where 

= £ + 

Si 

and let 

where the product is over all inversions ((3, a) in a permutation a G S>at and Sp a 
is the 5-matrix of the corresponding model. Also, we assume that a differential 
incorporates 

Theorem 3.1. The transition probabilities of the two-sided PushASEP, the ASAP 
and the AZRP are given in the form of 

P Y (X;t)=J2 I ■ ! ^IlSir" 165 * 1 '* 1 '"^ (3 ' 28) 

o-e§jv Cr Cr i 

where C r is a circle centered at zero with sufficiently small radius r so that all the 
poles of the integrand except at the origin lie outside C r . 

The proof of the theorem consists of showing that (13.281) satisfies (i) the forward equa- 
tion (12.261) for all X G Z , (ii) the corresponding boundary conditions (I2.14l) . (l2.23p 
or (I2.25p . and (Hi) the initial condition ( I2.27P in the physical region of the corre- 
sponding model. The proofs of (i) and (ii) are the same as the proof of the ASEP 
case [21]. For (Hi) it suffices to show that 

£/•••/ ^nC7' t "" 1 *--^ = °- 

a^id JCr JCr i 

In case of the ASAP, the ^-matrix of the ASAP (I2.24p is in the same form as the 
S'-matrix of the ASEP up to constants, hence the proof is the same as the proof of 
the ASEP case [24J. The 5-matrix of the AZRP has an additional factor besides 
the S-matrix of the ASEP and the 5-matrix of the two-sided PushASEP is in the 
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reciprocal form of the S'-matrix of the AZRP. Hence it is natural to use the techniques 
of the proofs of relevant lemmas in [24J in order to prove (in). Let 

J Cr J Cr ^ 

Now we present some lemmas corresponding to Lemma 2.1 - Lemma 2.4 in [23] . 

Lemma 3.2. Suppose that (3 appears in position (3 + 1 in a and the entries following 
(3 are greater than (3. Then 1(a) = for the two-sided PushASEP. 



Proof. With the assumption the permutation a has a unique inversion of the form 
(a, (3). Let 7 be the position of a. There are possibly other inversions of the form 
(a, 5) but no other types of inversions. The number of inversions of the form (a, 5) is 
a — 7 including (a, (3). In this case, the variables £ a and £p in the integrand appear 

as 

(a,6) 

where the product is over all inversions of the form (a, 5). Substituting 

& = rrV' (3 - 29) 

in the integrand and in the differential, the variable £ Q appears as 

cXj-xp+x+yp-ya+a-y 

and the exponent of £ Q satisfies 

a; 7 - xp+i +yp-y a + cx-l < -1 (3.30) 
because xp + i — x^ > ((3 + 1) — 7 and y a — yp > a — (3. Also we have 

TT at _ ( 1 + ^ IL^a,/? £,$ 1 ~ V Tl&^ip 1 

(11) a5 Ai+A^iL^e^-e, 



X 



n 



by the substitution f l3.29[) . If we integrate with respect to £ Q over C r , the order of 
the pole in 77 4T 1 * s g rea ter than 2 because of (I3.30p . Moreover, the singularities 
in the denominators are bounded away from zero because the second summands of 
denominators are 0(r 2 ) and \i 7^ 0, and thus 1(a) = 0. □ 
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The next lemma shows that Lemma 2.1 in [21] holds for the AZRP in the physical 
region of the AZRP. 

Lemma 3.3. (The AZRP version of Lemma 2.1 in [2$ ) Suppose that j3 appears in 
position — 1 in o and the entries preceding f3 are less than (3. Then 1(a) = for 
the AZRP. 



Proof. With the assumption the permutation a has a unique inversion of the form 
((3, a). Let 7 be the position of a. There are possibly other inversions of the form 
(5, a) but no other types of inversions. The number of inversions of the form (5, a) is 
7 — a including a). In this case, the variables £ a and £p in the integrand appear 

as 

(*,«) 

where the product is over all inversions of the form (5, a). Substituting 



in the integrand and in the differential, the variable £ Q appears as 



£X~ t -xp_ 1 +yp-y a +-y-a 

and the exponent of £ a satisfies 

x 7 - xp-x +y P --y a + r y-oc>l 

because x 7 — xp^i > 0, yp — y a > and 7 — a > 1. The rest of the proof is the same 
as the proof of Lemma 2.1. in [21]. □ 



Lemma 3.4. Let a G §at and (/3,a) be an inversion of a G S^r. Then, 

(fa) ^ < 



Proof. There are N — i entries to the right of a(i) and i — 1 entries to the left of o~(i) 
and o~(i) — 1 entries are smaller than cr(i). Suppose that there are n inversions of the 
form (o~(i),cr(j)). This implies that the number of entries smaller than a(i) is n and 
the n entries lie to the right of o~(i). Hence o~(i) — 1 — n entries less than o~(i) lie to 
the left of o~(i), which implies that (i — 1) — (o~(i) — 1 — n) entries are greater than 
a{i) and lie to the left of cr(i). □ 
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Lemma 3.5. Let x^yi G Z and (b,a) be an inversion of a permutation a G Sjy. 
Define 

and denote by the function obtained by substituting £ 7 = ^ J | in the function h a . 

Assume that Xi < x i+ i and yi < y i+ i for all i. (i) If two permutations a and a' differ 
only by an interchange of two adjacent entries a and (3, then h1\^ a= ^ p = h1,\^ a= ^ 
when 7 7^ a, (3. (ii) Moreover, if 7 is to the right of a and (3 and 7 < a, (3, then the 
largest power of £ Q and £p in h? a and h^, is zero. 



Proof. Let a < j3. Assuming that a(i) = a, a(i + 1) = j3, a'(i + 1) = a, a'(i) = (3 
and cr(j) = cr'(j) = 7. Then by Lemma 3.4 the variables £ a , £p and £ 7 appear as 

C%i— Vu — 1+a— i £ x i+i— Vp — £Xj— y 7 — 1+7— j 
Sa S/3 S 7 

in h a and appear as 

{Xi+i— y a — 1+a— c x i— ^ _1 +/3-Vi 3 - y 7 — 1+7-j 
Sa S/3 S7 

in hfjf . By the substitution (7 7^ «,/?), the variables £ a and ^ appear as 

tx i -x 3 +y 7 -y C( +a-7+j-jt :I: i+i~ :c j+?/7-^+/ 3 -7+i-(*+l) 
So S/3 

in ^ an d 

£Xi+i— Xj+y-y— y a +a— j+j — c Xi ~ x j Jr V'i~ yp+P~~l+3~ ' 

Set S/3 

in Thus, (i) follows from the fact that the sum of exponents of £ a and £p in /i^ is 
equal to the sum of exponents of £ Q and £p in hi,. Now, suppose that j > i + 1 and 
7 < a, (3. Then (ii) follows from the fact that Xi < x i+i < Xj and y 7 < y a < yp. □ 



The next lemma is for the AZRP. The proof is almost the same as the proof of Lemma 
3.5 and so we omit the proof. 

Lemma 3.6. Let Xi,yi G Z and (b,a) be an inversion of a permutation a G §/v. 
Define 

m*i,-,6t) :=nlnar w_1 

(M i 

and denote by W a the function obtained by substituting £ 7 = ^ - | in the function h a . 
Assume that x t < x i+ i and y^ < y i+ \ for all i. (i) If two permutations a and a' differ 
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only by an interchange of two adjacent entries a and (3, then h1\^ a= ^ B = h1,\^ a= ^ 
when 7 ^ a, (3. (ii) Moreover, if 7 is to the left a and (3 and 7 > a, (3, then the 
powers of £ a and £3 in hi and h^, are greater than or equal to 2. 

Lemma 3.7. Suppose that in the permutations a and a' the entry 7 appears to the 
right of two adjacent entries a, /3 > 7, and that the permutations differ only by an 
interchange of a and (3. Then 1(a) + I(a') = for the two-sided PushASEP. 

Proof. We will show that 1(a) + I(a') = when we integrate with respect to £ Q and 
£p. Suppose that a < (3, a = a(i), (3 = a(i + 1) and 7 = a(j) with j > i + 1. From 
the assumption there are inversions (0,7) and (/3,7) in both a and a' and there is 
additionally ((3, a) in a'. For a the variables £ Q ,£/3 and £ 7 appear as 

er^~ (4+1) er j ^ A II sUr^er 1 "^' 1 ^ 7 " 1 ( 3 - 31 ) 

in the integrand where the product is over all the other inversions involving a, [3 and 
7 excluding (a, 7) and (/?, 7). All the ^-matrices in a appear in a', too. Substituting 
£ 7 = j-^J ^ in the integrand and in the differential, the powers of and £g are less 

than or equal to —1 from Lemma 3.5. The product S^S^ is 

» + x v Us^n c 1 - fj- + a?? n«5^, 7 - & 

which is in the same form as the one in Lemma 2.2 in |24j. Thus, with the same 
argument as in Lemma 2.2 in |24j we have £ a = £p. Similarly, for a' the variables 
£ a ,£/3 and £ 7 appear as 

£Ot-(i+l)j;P-ijn-j ot ct Ct TT S t x !~y/3~ 1 f3;i+i-i/c ( -lt2,j-y 7 -l (o oo\ 

Sa ?/3 S 7 °/3a°a 7 °^ 7 11 °<5<5"?/3 Sa S 7 l^O.OZj 

in the integrand. Substituting £ 7 = j-^J ^ , it follows from Lemma 3.5 that ( 13 .3 If) is 

the negative of (I3.32p when £ a = because Sl a = — 1 when £ a = £p. Hence the 
integrands in J(cr) and /(cx') differ only by the sign after integrating with respect to 
£ a and £0 and so /(a) + I(a') = 0. □ 

The following lemma can be proved by using Lemma 3.6 and the same procedure as 
in the proof of Lemma 3.7, and hence we omit the proof. 
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Lemma 3.8. (AZRP version of Lemma 2.2 in [24\j ) Suppose that in the permutations 
a and o' the entry 7 appears to the left of two adjacent entries a, (3 < 7, and that the 
permutations differ only by an interchange of a and (5. Then 1(a) + I(a') = for the 
AZRP. 

While Lemma 2.4 in [21] is for the ASEP and the AZRP, the next lemma is for the 
two-sided PushASEP 

Lemma 3.9. (The two-sided PushASEP version of Lemma 2.4 in WM) For any N 
the set E>n \ {id} is the union of disjoint subsets, each of which consists of either of a 
single permutation satisfying Lemma 3.2 or a pair of permutations satisfying Lemma 
3.7. 

Proof. It is clear that the result holds for N = 2. Assume that the result holds for 
N—l. Let A a be the set of permutations in which a appears in position N. If a = N, 
we apply the induction. Notice that all permutations in A^-i satisfy Lemma 3.2. If 
a < N — 1, we can pair (N — 1)! permutations in A a as in Lemma 2.3 in [24J, and 
these pairs satisfy Lemma 3.7. □ 

Hence, for the two-sided PushASEP the proof of (in) of Theorem 3.1 follows from 
Lemma 3.2, 3.7, and 3.9 and for the AZRP the proof follows from Lemma 3.3, 3.8 
and Lemma 2.4 in [21]. This completes the proof of Theorem 3.1. 

An immediate result of Theorem 3.1 is a natural mapping between the state spaces 
of the ASEP and the AZRP under which transition probabilities are conserved. This 
can be considered as another approach of the well-known mapping between the ASEP 
and the AZRP [6l [9l [12j [21] . Specifically, let S and S' be countable state spaces of the 
ASEP and the AZRP with N particles, respectively, and X, Y e S and X', Y' e S'. 
Let A^ and A^ be coefficients of Bethe ansatz solutions of the ASEP and the AZRP, 
respectively. Then under the bijection / from 5" to S defined by 

/ : (Xi)i ->■ (Xi + i)i, (3.33) 
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we have from Lemma 3.4 



ctGSjv 

E 
E 

Py(X;t). 



Cr 



(i) 6 



N 



x{+i— y' cr(i)— 1 \j 



Cr 



A a II 

i 



N 



Hence the distributions of the mth particle's position from the leftmost in the ASEP 
and in the AZRP are conserved under the mapping (I3.33j) because the distributions 
of the mth particle's position are obtained by summing transition probabilities over 
all possible configurations. Thus, F(x m (t) = x), the probability of finding the mth 
particle at x at time t in the AZRP is given by (5.12) in [23] with t/i replaced by yi + i. 
Let us modify the iV-dimensional integrand (1.1) in [23] by letting 

£j — Ci 1 — £i • • • £;v 



I z (x,Y,S) = H 



i<j 



v + i^j - Ci (i - 6 



N, 



n e 



-x-( yi +i)-l B fc) t 



and as in [23], denote by Iz(x,Y$,£) the integrand analogous to Iz(x,Y,£), where 
only the variables £j with i G S occur. Then we can immediately write the explicit 
form of P(x m (t) = x) for the AZRP. 

Proposition 3.10. Let S C {1, • • • , N} and a(S) = Y. ie s *■ Then 

F(x m (t)=x) = (-l) m +\p q ) m ( m -W (3.34) 



|S|>m 



\S\ -1 

l^l — m 



P 



a(S)-m\S\ 



Q 



<r(S)-|S|(|S|+l)/2 



Cr 



I z (x,Y s ,0d ]s ^ 



Cr 



where Cr is a circle centered at with sufficiently large radius R so that the poles of 
the integrand are in the circle. 



It is natural to expect that ¥(x m (t) = x) in the ASEP with step initial condition 
should be equal to that in the AZRP with the delta function-like initial condition 
that infinitely many particles occupy the origin and all other sites are empty. This 
can be simply verified as follows. Let S = {zi, ■ ■ ■ ,Zk} in (13.341) . As in the proof of 
Corollary in [24], by the change of variables, the factor 
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in the integrand Iz(x,Yg,C) becomes 



i 

for the delta function-like initial condition. Hence the rest of the proof is the same 
as in the proof of Corollary in [21] . Moreover, the AZRP with the delta function-like 
initial condition has the same asymptotic behavior as the ASEP with step initial 
condition as well as the Fredholm determinant representation [251 12S] . 

Remark 1. The S-matrices of the ASAP and the ASEP are equivalent up to con- 
stants. Specifically, if p and q in the S-matrix of the ASEP are replaced by —///A and 
1/A, respectively, we obtain the S'-matrix of the ASAP. Hence the integral formulas 
for the ASAP are obtained from the integral formulas of the ASEP with p = — fi/X 
and q—l/X and without any change in the energy term, e(£i). 

Remark 2. While there is a direct mapping between the ASEP and the AZRP, on 
the other hand, there is no such mapping between the ASEP and the two-sided 
PushASEP. Moreover, there is no such simple relation between the S'-matrix of the 
PushASEP and the S-matrix of the ASEP as the ralation between the S-matrix of 
the ASAP and the S-matrix of the ASEP. Thus, it needs further investigation to 
study the current distribution and its asymptotics in the two-sided PushASEP and 
they remain as problems for the future. 
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